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We show that controllable inhomogeneous coupling between two-level systems and a common data bus pro- 
vides a fast mechanism to produce multipartite entanglement. Our proposal combines resonant interactions and 
engineering of coupling strengths — between the qubits and the single mode — leading to well defined entan- 
gled states. Furthermore, we show that, if the two-level systems interact dispersively with the quantized mode, 
engineering of coupling strengths allows the controlled access of the symmetric Hilbert space of qubits. 



Introduction. — Nowadays, complex quantum architectures 
are being developed aiming at physical implementations in 
quantum information processing fP-H^. Their complexity is 
often revealed by intrinsic inhomogeneities which may limit 
the quantum control, for example, to prepare entangled states. 
This situation is found in a quantum dot, where the spatial 
distribution of the single-electron wave function makes its 
coupling to the spin bath position-dependent f7|. Inhomo- 
geneities also appear in circuit QED fSllSl, where the cou- 
pling strength of light-matter interaction depends on intrinsic 
parameters of qubits |6, 7 |. Other physical scenarios exhibit- 
ing inhomogeneities are impurities embedded in silicon O, 
as well as hybrid systems involving NV-centers O [TOl . 
Nonetheless, the advent of these platforms — particularly cir- 
cuit QED — brings new possibilities for engineering direct or 
indirect coupling among qubits in order to produce multipar- 
tite entangled states 1 11|. Most implementations are based on 
slow off-resonant interaction among the qubits and a common 
data bus, however, resonant interactions have the advantage 
that their gating time scales with the inverse of the coupling 
strength, resulting in faster two-qubit gates |[T2ll . 

In this work, we show that engineering the inhomogeneities 
of the coupling strengths between the two-level systems and a 
single-mode resonator, provides a fast mechanism for gener- 
ating multipartite entangled states. First, we discuss the sim- 
plest case of a two two-level system interacting inhomoge- 
neously with a single-mode resonator, looking for the specific 
conditions to generate well defined amounts of entanglement. 
We extend this study to the multipartite case, where our pro- 
posal allows the generation of W entangled states. Second, 
we consider a dispersive interaction among qubits and the 
single-mode resonator such that we the engineering of cou- 
pling strengths allows the controlled access of the symmetric 
space of qubits. In particular, we profit from the resulting 
quantum dynamics to propose a mechanism to produce differ- 
ent classes of entangled states. Finally, we present our con- 
cluding remarks. 

Two-qubit system interacting inhomogeneously with a sin- 
gle quantized mode. — Let us start our discussion considering 
the simplest case of a two-qubit system interacting with a sin- 
gle quantized mode. A main assumption of our model will be 
the individual addressing of qubits, allowing the controllabil- 
ity of their coupling strengths with the single mode. We aim 



to generate correlated states between qubits, so the physical 
mechanism has to be the exchange of a single excitation me- 
diated by the quantum mode. If the first qubit is in the excited 
state, the second qubit in the ground state, and the quantized 
mode having no excitations, the single excitation located in 
the first qubit will be transferred among the parties according 
to the Hamiltonian 

H = hxi{\ei){gi\d^\gi){ei\d'^)^hx2{\e2){g2\a^\g2){e2\a'^), 

(1) 

where each qubit is coupled inhomogeneously (xi and X2) 
to the quantum mode. As the initial condition is lipo) = 
|ei5'2)|0), the accessible Hilbert space for the composed sys- 
tem is spanned by vectors {|ei^2)|0), |^i^2)|l), |^ie2)|0)}, so 
that the state of the system reads 

iV^(t)) =Ci(t)|ei^2)|0)+ C2(t)|^1^2)|l)+C3(t)|^l 62)10). 

(2) 

The quantum dynamics supported by Hamiltonian ([T]) in- 
duces entanglement between the qubits through their inter- 
action with the quantum bus. However, it would be notice- 
able to produce maximal entanglement between the qubits in 
a single step. From Eq. (|2]), a maximal entanglement could be 
achieved when the probability amplitude C2{t) vanishes. The 
probability amplitudes in Eq. ^ are given by 

ci(t) = l + ^[cos/it-l] (3) 

C2(^) = —i^ sin /it (4) 

C3 (t) = [cos /it - 1] , (5) 



where /i = vXi+xf • Imposing the condition C2 (t) = 0, we 
obtain ci{t) = 1 - 2x?//i^ and cs{t) = — 2x1X2/^^- Thus 
resulting in the entangled state 

m = tt/m)) = (1 - 24)|ei52> - 2^15162). (6) 

/i^ /i^ 

We realize that a maximally entangled state cannot be 
achieved in a single step with resonant interactions, and both 
qubits having same the coupling strengths. Here the inhomo- 
geneous coupling will play a key role. In order to reach a 
maximally entangled state \ip{0)) = {\e1g2) + \gie2))/V^, 



2 



the coupling strengths, xi and X2, have to satisfy the relation 
X2 = -(V^+l)xi orx2 = {V2-l)xi- This simple analysis 
shows that a suitable control of couplings strengths can be a 
useful mechanism to control the access to the Hilbert space of 
a two-qubit system. The natural extension of this idea is the 
case of a A^-qubit system. 

Consider an array of N qubits, interacting in- 
homogeneously with the quantum mode. At resonance, 
this situation can be described by the Hamiltonian 



H 



\g^){e^\a^)• 



(7) 



Let us sssume the k-th qubit initially in the excited state and 
the remaining qubits be in the ground state, that is, \ip{0)) = 
\gg..ek--gg). For this initial condition, the evolution given by 
Hamiltonian ^ leads to the state 



=co(t)|^i^2..)|l) + ^Q(t)|^i^2..ei^i+i..)|0), 



(8) 



where the probability amplitudes read 

■ Xk 



cq = —I — sin fit^ 

Cj - djk + 



COS /it — 1], 



(9) 



with ji = Y Xp Cfe(O) = 1 and Cj^/c(0) = 0. From these 

expressions we learn that inhomogeneous coupling allow us 
to achieve a maximally multipartite entangled 1^- state reso- 
nantly and in a single step, by setting the interaction time to 
t = it/ ji and 5jk — 2xjX/e/M^ = 1/vW. This leads to the 
following condition for the coupling strengths 



XI j Xj 



(1 - ^N) 



1 



(10) 



For example, if the first qubit is initially excited and has a 
coupling constant Xi, and the remaining qubits being in the 



ground state with a coupling constant x, the condition ( 10) is 
fulfilled when xi = (1 + V^)X' resulting in the W state 



= 7r//i)) = --=^\gg ...Cj ...gg)\{)). 



(11) 



Thus, we have been able to find conditions under which in- 
homogeneity allow the preparation of maximally multipartite 
entanglement with a single resonant interaction. Nonetheless, 
this physical mechanism is limited to the generation of en- 
tangled states involving a single excitation. Naturally, other 
limitations are those imposed by dissipative effects acting on 
the qubit or the quantum mode. 

A physical setup where these ideas could be implemented, 
considering the available experimental resources, is circuit 
QED ISHll. Two-qubit interactions mediated by a single 
mode of a resonator have been reported by several experi- 
mental groups |[T3UT6ll . In particular, in Ref. 1 14], the Tavis- 
Cummings model 1 171 has been demonstrated experimentally 
up to three qubits, reaching maximally entangled states. 
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FIG. 1 . Fidelity of maximally entanglement as a function of the num- 
ber N of qubits. We have considered that the first qubit has the 
experimental parameters following Ref. 1 18|, while the remaining 
ones follow the corresponding engineered coupling values. Maxi- 
mum coupling: xi = 54 MHz, cavity k = 3.2 MHz and the qubit 
spontaneous emission rate 7 = 0.6 MHz. The initial condition is 

|V^(o)) = |^^...^)|i). 



As an example, in Fig. [T] the fidelity of the preparation 
of a TV-particle W state in a circuit QED setup is shown as 
a function of the number of qubits. In this simulation, we 
have considered dissipative mechanisms for cavity and qubits 
(dephasing and relaxation) using parameters of recent experi- 
ments 1 18 1. It is clear that the fidelity of the W state increases 
as a function of the number of qubits, as expected from the 
higher effective coupling ~ \/N x- 

Another physical setup where maximal entanglement can 
be reached resonantly in a single step is trapped-ion sys- 
tems |[T9l . Here, an alternative procedure can be followed: 
consider the first qubit tuned in the blue sideband transition 
and second qubit tuned in the red sideband one. For an initial 
condition |V^(0)) = \g1g2), the system will evolve within the 
subspace spanned by {|^i^2)|0), |ei^2)|l), |eie2)|0)}. Im- 
posing similar conditions as for the previous case, this pro- 
cedure will lead to the maximally entangled state \ilj{t = 
/i^)) = {\e1e2) + \gig2))/V2). 

Collective Control of Symmetric Dicke space. — A resonant 
interaction between qubits and the quantum mode limits the 
possibility of generating other states than a W state among 
qubits. As the quantum mode acts as an intermediary system, 
we can explore a regime where it can be adiabatically elimi- 
nated, allowing a direct transfer of excitations between qubits. 
In this way, a set of qubits having a given number of k excita- 
tions could directly interact with another set of qubits having 
q excitations, given rise to the generation of states with k -\- q 
excitations. To this purpose, we consider a collection of N 
qubits, coupled to a quantum mode with coupling strengths 
gi and detuned from the qubit transition by Ai. In addition, 
we consider M qubits coupled to the same quantum mode 
with coupling strengths ^2 and detuning A2. This setup is 
schematically shown in Fig |2] The Hamiltonian describing 
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FIG. 2. (Color online) Schematic representation of N qubits coupled 
to a single mode with coupling strength gi and other M qubits with 
a coupling strength g2. 



this situation in the interaction picture reads 



H.C., 



where we have defined the collective operators 



(12) 

Ej^i |e)j(g| and Jt = |e)j(g|. In the far off- 

resonance regime, i.e. Ai, A2 ^ ^1, ^2, an effective interac- 
tion between the A^-qubit and the M-qubit subsystems arises. 
Under the assumption that no photons are initially present in 
the quantum mode, the effective Hamiltonian reads 



^eff = nXiS^ s ^ n\2J^ J 



where A,- = g'^J/S.j, (^eff = ^1^2/ A and A 



'0 - I 
A2). By noticing that 

N 



2AiA2/(Ai- 



N 



S^S= aja,+^|e,.){e,.|, (14) 

M M 

the Hamiltonian ( p3]) can be rewritten in a more convenient 
form /^eff = ^0 + Hj, where 

N M 

Ho = ^Ai ^ + hX2 ^ 



Hi — hQef[ 



(16) 



where J = A2 — Ai + A2 — Ai is the effective detuning pa- 
rameter. This Hamiltonian is the starting point for the sub- 
sequent discussion about the selective control of symmetric 
Dicke space. Indeed by conveniently choosing initial condi- 
tions for the set of N and M qubits we can force the evolution 
within symmetric Dicke space. In that follows, we will con- 
sider the general problem of the evolution of the set of N and 
M qubits, prepared in a special product of symmetric Dicke 
states. The same protocol will show us how to prepare a set of 



qubits in such states. Let us consider the overall system ini- 
tially in a product state of the A/'-qubit subsystem in a symmet- 
ric state with k excitations and the M-qubit subsystem having 
q excitations, that is 



1*0) 



IDDlDf), 



(17) 



where |D^) stands for a A''-particle Dicke state with j exci- 
tations. In general, this initial state will evolve — under the 
action of the Hamiltonian — within the subspace 



{|Dn|Df),|D^±i)|Df^i),|D^±2>|Df^2),---}- (18) 

We aim to gain control of the symmetric space by reduc- 
ing the size of the effective Hilbert space that the system can 



visit. More precisely, given the initial state ( 17 ) we would like 
the evolution be restricted to the two-dimensional subspace: 
{|Df )|Df ), |Df±i)|D,%i)}. This can be done by adjusting 
the parameter S, setting this transition to resonance and leav- 
ing other possible transitions far from resonance. To see how 



this works let us analyze how the Hamiltonian ( 16 ) acts on the 
initial state ^Tl) 



with 



Hoi's 0) 



fk,q 



^^eff(//c+l,g-ie" 



ID 



M\ 



ID 



M 



Xik{N -k)■ 



l){N-k){M 
-\2q{M -q). 



(19) 
(20) 



Thus, each transition {|Df )|Df ), |Df^i)|D^i)} has asso- 
ciated a detuning equal to S = — {^k±i,qT'i^ — ^k,q) • Setting 
(5 = for given values of k, q, N and M, for a chosen transi- 
tion other subspaces will be far from resonance, provided that 
this detuning is much larger than the effective coupling (^eff- 
As an example, let us consider the case when the first subsys- 
tem has single qubit in the excited state (A^ = 1 and k = 1) 
while the second subsystem is initially in a symmetric state 
with M qubits and q excitations, that is 



M 



|^o) = |e)^a,|Df). 



(21) 



Considering this initial state and setting the parameters such 
that 6 = S — Si^q-^i — 5o,g = 0. we enforce the system to 
evolves only in an effective two-dimensional Hilbert space 
spanned by {|e)|D^),|^)|D^^)}. This selectivity is pro- 
vided by S r^eff, such that others subspaces do not evolve. 



When these requirements are fulfilled, the initial state ( 22 ) be- 
comes 



1^. 



e) |Df ) + i sin (l^eff /o,,+it) Iff) |D^i)] . 



For example, in Fig.|3]we show the dynamics for the case 
M = 6 and q = 3. Here we see that the dynamics occurs only 
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FIG. 3. (Color online) Evolution of populations of states composing 
the initial condition: |^o) = (1/^)1 e) ^f^^ aj\Df) for M=6 as a 
function the time r = 2/o,4^^efft/7r. The initial probability ampli- 
tudes are ao — 0.1, ai — 0.2, a2 — 0.3, as — 0.4, — 0.5 and 
6 is a normalization constant. The parameters are: gi/g2 = 0.02, 
A2/A1 = 1.002 and Ai/^2 = 100. 



in the selected subspace spanned by {|e)|D3^), |^)|D4^)}, 
while the population of other states remaining virtually un- 
changed. This reduced dynamics can be used to prepare arbi- 
trary Dicke states by measuring the state of the first qubit |20|. 
That is, if we allow the system to evolve during a time t = 
7r/(2(^eff/o,g+i) and measuring the state of the first qubit we 
will find that the M-qubit subsystem will collapse into a Dicke 
state |D^i) with probability la^p. 

However, preparation of arbitrary Dicke states can also be 
achieved deterministically using the ideas above. For exam- 
ple, we consider the system initially in the state 



1*0) = |e>|D, 



(22) 



Under the action of Hamiltonian ( [161 ), this initial state will 
couple such as 

^o|^o) = Mi,,|e)|Df) 

leading to an effective two-level dynamics in the subspace 
{|e)|D^), |^)|D^i)}. If we set the parameter such that, 
^ = — ^o,g+i. the evolution of the initial state will be 



that is, it becomes in a (M + 1) -particle Dicke state with q-\-l 
excitations. 

Following this idea, arbitrary Dicke states can be deter- 
ministically prepared in a sequential manner. This works 
as follows: Consider and initial state of the form |^o) = 
|e)|Do^) = \e)\gig2 . . . ^at). Then, if the system evolves dur- 
ing a time ti = tm,o^ we will have l^t^) = \D^~^^). Now, 
we consider an additional qubit initially in the excited state, 
that is, l^^ti) |e)|Df^+^). If this system evolves under the 
action of Hamiltonian (16) during a time t2 = tm+ 1,1, the 
system will be in the state l^ts) = In this manner, 

following this procedure we can sequentially prepare after q 
steps, the Dicke state |D^+^). 

It is worth to mention that different proposals to prepare 
Dicke states have been presented for homogeneously coupled 
systems ll2QH24l and in inhomogeneous systems in a non- 
deterministic procedure ll25ll . 

Entangled Dicke States. — In the following we will show 
that this quantum dynamics can be used to prepare a different 
class of entangled Dicke states. Let us consider the situation 
described in Fig.|2]for N = M. In such case we look for the 
conditions required to achieve an effective two-level dynam- 
ics in a subspace {|Df )|D^), |Df^i)|D^^i)}. It is clear that 
under these conditions we should be able to generate entangle- 
ment between Dicke states. A particular state we are able to 
prepare using the procedures described above is a NOON-like 
state of the form 



|*) = -i.(|D]^)|Do^) + |Do^)|D]^)). 



(26) 



This can be done as follows. First, we consider the ini- 
tial state 1^0) = |D]v)I^o^)- Then, we set into resonance 
the transition in subspace {|D]^)|D^), |D]^_i)|Df )}. Set- 
ting the interaction time such as t = 7r/(4(^eff/Ar-i,i), the 
state will evolve to 



1*1) = 



1 

71' 



IDiy\iT^iv\ 



i)|Df))- (27) 



By setting parameters such 
{\T)%_,)m,\D%_,)\T)^)} is 



that the subspace 
J ^±^-^ /, I j_y^_2 / 1 -i->'2 // resonant, the state 

|D]^)|Dq^) does not evolve since its associated transi- 
tion is far from resonance. We lead the system to evolve 
during a time t = 7r/(21]eff/Ar-2,2), obtaining 



1^,) = cos (^]eff/o,g+it) |e) \D'^')^i sin (l^eff/o,,+it) \g)KU) 

When this state evolves during a time t = 7r/(2(^eff/o,(?+i). 
the resulting state will be a symmetric M-particle Dicke state 
with q-\-l excitations 



1*2) = 



^i\I)^)\D^) + e\I)%_,) 



D^)). (28) 



Along the evolution we transfer the excitations from the first 
set of qubits to the second set of qubits in such a way 



= \9)\K+i)- 
Now, if the interaction time is set to be 



1 



^eff/o,g+l 

the system will evolve into 



arcsm 



M-q 
M + 1' 



(23) 



(24) 



(25) 



|D]^-2>|D2 



N\ 



ID 



N 

N-3 



Do^)|D-.) (29) 



In this manner, after N steps the state ( |26| ) can be prepared. 
As an example, we consider the case of each system having 
N = 3 qubits. The required values for the parameters to set 
into resonance the needed transitions — to prepare an entan- 



gled state of the form ( [26| ) — are listed in Table I. 

The scheme for finding the exact values to select a given 
transition is as follows: first, as pointed out, we must fulfill the 
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Step 


6 


A2/A1 


91/92 


^1/91 


^2/92 


I 


3(Ai + A2) 


0.9996 


18.630 


53.7 


1000 


II 





1.0025 


70.621 


20 


1416 


III 


-3(Ai + A2) 


1.0075 


17.611 


20.15 


355 



TABLE I. Parameters for the preparation of the state ([26| 
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FIG. 4. (Color online ) Schematic of sequential preparation of 
NOON-like state in Eq. (pOl. 



following requirement, Ai, A2 ^ ^1, ^2- To tune the selected 
transition we must set the parameters gi, Ai, A2, ^2 or ^2 in 
such way that the relation J = is satisfied. To achieve this, 
we may choose either gi as a function of Ai, A2,^2 or ^2 
as a function of Ai, A2,^i. The choice of these parameter 
must be such that 6 (^eff for the other possible transitions. 



preventing in way, evolution in subspaces different from the 
chosen one. A direct consequence of such requirements is that 
the qubits cannot be homogeneously coupled, that is, gi ^ g2- 



In Fig.|4] we show the process to prepare a multipartite en- 
tangled state 1^) = (l/V2)(|Df )|D^) + |D^)|Df )) for 
= 3. The first step — as shown in Fig. |4] — is to prepare 
a superposition 



|*) = (l/v^)(|D3^)|Dn + |D^>|Dn) 
The following two steps allow the transitions 



(30) 



|D^)|Df ) |Df )|D^) ^^5^ |D^) 



D 



N\ 



without changing the population of the state \Dq^ ) |Dq^ ) . 

Conclusion. — We have described the role of the inhomo- 
geneous coupling between a single quantized mode and a 
set of two-level systems. In particular, we have shown that 
combining resonant interactions and engineering of coupling 
strengths allows a fast mechanism to produce multipartite en- 
tanglement. Furthermore, in the dispersive regime of light- 
matter interaction, the versatility in the engineering of cou- 
pling strengths provides a mechanism for controlling the ac- 
cess to the symmetric Hilbert space of qubits. At the same 
time, it allows the access to special classes of entangled states, 
both probabilistically and deterministically. 
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